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Real and Image Fields of a Relativistic Bunch 
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Skobeltsyn Institute of Nuclear Physics, Moscow State University 

Abstract 

^O I We derive analytical expressions for external fields of a charged relativistic bunch 

^^ ' with a circular cross section. At distances far from the bunch, the field reduces 

^SJ ■ to the relativistic modified Coulomb form and in the near region, reproduce the 

Vh ! external fields of a continuous beam. If the bunch is surrounded by conducting 

^^' surfaces, the bunch self-fields are modified. Image fields generated by a bunch 

"^ ■ between two parallel conducting planes are studied in detail. Exact summation 

C^ I of image fields by the direct method invented by Laslett allows the infinite series 

,__i' to be represented in terms of elementary trigonometric functions. 

^: 
6 ■ 

c^ ■ 1 Introduction 

O ■ In an accelerator, the charged beam is influenced by an environment (beam pipe, accelerator 

c/3 , gaps, magnets, collimators, etc.), and a high-intensity bunch induces surface charges or 

^^1 currents into this environment. This modifies the electric and magnetic fields around the 

Q_^' bunch. There is a relatively simple method to account for the effect of the environment by 

'~~', introducing image charges and currents. 

^— I I Over forty years ago Laslett PP analyzed the influence of the transverse space-charge 

K^ ' phenomena, due to image forces, on the instability of the coherent transverse motion of an 

I . intense beam. Methods of image fields summation are described in his paper ^ which pre- 

f^ I sented some field coefficients calculated for infinite parallel plate vacuum chambers, magnetic 

^r ■ poles and vacuum chambers with elliptical cross sections and variable aspect rations. The 

vQ ■ resulting image field was calculated only in the linear approximation and depends linearly on 

^^ i the deviations x and x of the bunch center and the position of a test particle, respectively, 

lyj '_ from the axis. They act therefore like a quadrupole causing a coherent tune shift. The 
approximation used is incorrect if the field observation point x is located far from the bunch 



CJ 



C^ 



^, or if the bunch center is close to a conducting wall. 

r~| I In the present paper we consider the problem of the image field summation once again for 

Oh' a very simple geometry, namely, a relativistic bunch moving between infinitely wide parallel 

conducting planes. The problem is far from being a pure academic one J2j. In applications, 
k><( I in particular by study of the electron cloud effect |H] and the dynamics of photoelectrons in 

^ ' the beam transport system, it is important to know the distribution of electromagnetic fields 

not only in vicinity of the bunch, but in the whole gap. We did not find publications with 

attempts to suin up the series (|17|) in an approximation beyond the linear one. In Section 3 

we present the exact solution of the problem. 

To solve the problem formulated above, in Section 2 we first derive an equation for the 

external electromagnetic field generated by a cylindrical bunch of charged particles. The 

task is specified as follows. 
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The external radial electric E±, and azimuthal magnet B^ self fields for a round un- 
bundled relativistic beam of the radius b and a uniform charge density are [3]-[7j 

E. = .'A (1) 

r 

Ho 2qX 

B4, = -. c/3, 2) 

Att r 

where n = l/47reo, A is the linear charge density, q is the charge, /5 = v/c is a normalized 
velocity of the beam constituents and c is the velocity of light. In many applications, equa- 
tions (^ and are used to describe fields of an individual bunch too. However, in the form 
(^, the bunch fields do not depend on the bunch energy and at large distances do not 
follow the Coulomb asymptotic. This contrasts sharply with the fields produced (at t = 0) 
by a rapidly moving single charge q 



E = K—^ 



1-/32 



1-P^sin^e 



3/2 f - - - 

-, B ^ (3x E, (3) 



where 9 is the angle which the vector r makes with the z-axis. Along the direction of motion 
the electric field is become weaker in 7^ times, while in the transverse direction the electric 
field is enhanced by the factor 7 

^^ = -7?- (4) 

Here, 7 denotes the particle Lorentz factor. 

In the next section we derive an expression for the transverse component of the bunch 
electric field, which the defects indicated above are rectified, and find the conditions at which 
the bunch fields are represented by ^ and ©. 

2 Self-Fields of a Charged Finite Cylinder 
with a Circular Cross Section 

Let us consider a bunch of charged particles uniformly distributed with a density p within a 
cylinder of length L and an elliptical cross section. The ellipsoid semi-axis in the x-y plane 
are a and b and the coordinate z-axis is along the bunch axis. Suppose that the bunch is 
moving along the z-axis with a relativistic velocity v = cf3. 

To compute the radial electric field of such a rapidly moving bunch, we have to sum up 
fields of the type ©, generated by the bunch constituents. In this way we get ji] 



E^{r, e, z) = Kp^{zh + {L- z)l2} (5) 



with 



J ^ f f {r - a cos{^ - (I))) adadcp 

J J (r^ + a^ — 2r(TC0s(^ — <f)))y^^^z^ + r^ + a"^ — 2racos{S, — <j)) 

J ^ f f {r - a cos{^ - (p)) adadcf) 

J J (r2 -hcr2 - 2ra cos{^ - (p))s/'y^{L - zy + r^ + a"^ - 2mcos(^ - </>) 

where a is the distance in the x-y plane from the z-axis to the elementary charged volume 
and 

< a < " < (/>< 27r. (8) 

y a^ sin (p + b"^ cos^ (j) 



Equation © represents the radial electric field as observed at a distance r from the bunch 
axis, at an angle S, relative to the x axis and at a distance z from the bunch tail. 

The integrals /i and I2 can be estimated only numerically 4^, if integrands are taken as 
it is. However, the integrands are easy to simplify if the bunch is relativistic, 7^1, and we 
would like calculate the field in vicinity of the bunch, r ^ L, but at distances much larger 
than the bunch radius, b <^ r. 

To simplify, we make use of the notation 

1 



r 



B = Acos{^-(l)), Y = A^-2B, Ci 
and the integrand of Ji can be written as 



1 + 



72^2 



X = Ci-Y 



(r- + YzY^^^M^ -B){1 + Y)-\l + X)^^/2 



..2 2\ 



(9) 



Now we expand the above expression in a power series by using A as a small parameter and 
keeping only terms up to the power A^ at each step. For the bunch shaped as a circular 
cylinder, a = b and we may set ^ = 0. Due to the fact that 

I-2TV 



cos 



-,2fc+l 



0, 



(10) 



all odd power of B vanish after integration in (j). This greatly simplifies the series generated 
from ®. After lengthy algebraic manipulations with Q, we get 

1 - (1 + lci)A^ + (2 + Ci + lcf)B^'^ 



(^2 + ^V)-i/2^ 



Substituting this expression in ©, we get 

7r62 



,62' 



By changing z^ to (L — z)^ in 1)121) . we obtain for I2 the following result 

7r62 



h 



where Co 



1 + 72(L - zf/r^ 



r\Jr'^ + 72(-L — z)2 
-1 



1 



S /)2 



(11) 



(12) 



(13) 



Notice that for particles uniformly distributed in the 

bunch volume, p = qN/irb'^L, where N is number particles per bunch. Substituting equations 
^-^ in ©, finally we arrive to 



E±{r,z) 



Lr 



1 + 



3 62 



ct] + 



1 + 



3 6' 



-r2 



)}■ <") 



V^r2 + 72^2 V" ■ 8r2-V ' ^r2+72(L 

This equation describe the electric field produced by a rapidly moving circular bunch. 

The field of a relativistic bunch described by H14|) , has different behavior at distances far 
apart of the bunch and in the near region, r < L. At very large distances, r ^ 7Z and 
r 3> 7(i — z), equation 1)14(1 reduces to the Coulomb form @. At the same time, in the near 
region and beyond the bunch tails, 72; ~ 7(L — z) ^ r and equation ()14|) simplifies to 

2QiVl 



E^ 



L r' 



which coincide with the external field (Q) of a continuous beam with A = N/L. 
Similarly we can show that the azimuthal magnetic field of the bunch is 

B, = --E^{r,z). 



(15) 



(16) 



3 Fields from Image Charges 

Following Laslett [1,^ (see also fSI), we consider a relativistic bunch of the length L between 
infinitely wide conducting planes at x = it/i. Suppose that constituents of the bunch are 
positively charged. For full generality, let the circular particle bunch be displaced in the 
horizontal plane by x from the midplane (0,y,z), and the observation point of the field be 
at (x, 0, 0) between the conducting parallel planes. The end points of the bunch are at 
z = zbL/2. The boundary condition for electric fields is Ez{ih) = on the conducting plane 
and is satisfied if the image charges change sign from image to image. Suppose that the 
distance between planes is of the order L. Thus, the electric field of each image is described 
by ((T3|) . To calculate the image electric field E^^i^nageix) in front of the plate, we add the 
contributions from all image fields in the infinite series ^ 

p . -. _ 2giV 

{2h-xiy^ -(2/i + xi)"^ -(4/i-X2)"^ +{4h + X2y^ 
+ {6h-xiy^ -{6h + xi)-^ -{8h-X2y^ +(8/i + X2)~^ 
+ (lO/i - xi)'^ - (lO/i + xi)-^ - {12h - X2r^ + (12/i + X2y^ + •••}, (17) 

where xi = x + x and X2 = x — x. These image fields must be added to the direct field of the 
bunch (jlSj) to meet the boundary condition that the electric field enters conducting surfaces 
perpendicularly. 

In the original paper pQ, the series ()17() was summed up only in the linear approximation 
in X and x, 

E±,image{x,x) = K——--^{2x + x). (18) 

The coefficient ei = 7r^/48 is known as the Laslett coefficient (or form factor) for infinite 
parallel plate vacuum chambers and magnetic poles. The approximation used in (|18)) is 
incorrect if the deviation of the bunch center from the axis is large (x ^ h) or if the field 
observation point x is located far off the bunch. Therefore, below we present the exact 
solution of the problem. 

In Appendixes A and B we prove that the exact summation of the series H17|) gives 

AqN 

E±,image{x,x) = K^-j^A{6,6), (19) 

where the image field structure function A depends only on normalized variables 6 = x/h, 
6 = x/h in the form 



TT cos{^5) 1 



A(M-)4[|- .,,7^'V^i -r4l- p°) 

2 12 sm(^()) — sm(^()) 0— oJ 

In Appendix A it is shown that in the linear approximation equation ()19|) recovers the part 
(fTSJ) derived by Laslett. 

We shall now estimate values of the function A in several particular points. If the obser- 
vation point of the field is located at the plane, x = h, then 6 = 1 and the structure function 
depends only on the bunch center position between planes, 6. Thus, from (|2()|) we get 



Equation (|21j) is singular at (5 — > 1 and shows that the conducting plane attracts the bunch 
with increasing force with the bunch displacement from the midplane. This phenomenon, 
involving the transverse movement of the bunch as a whole, arises from image forces and 
could lead to a transverse instability. 

For a bunch in the midplane, 6 = 0, the summed image field at the surface equals 






E±,^mage{h,0) = K-^A(1,0) = ^^^iri^ " l)' (^2) 



The image field l\W\i must be added to the direct field of the bunch ((T3|) to meet the 
boundary condition. Thus 



E±,tot{x, X) = E±^bunch + E±^image = l^-TT ( 1 + ^^(5, 5) ) . (23) 



2qN /I 
~Lh\5 



For a bunch in the midplane, 5 = 0, we find from 1)231) the expression of the transverse 
component of electric field generated by a relativistic bunch moving between wide conducting 

parallel planes 

2qN 7r/2 , , 

Lh sm(-|o) 

That is, at the surface, 5=1, the field is enhanced by a factor 7r/2 due to the presence of 
the conducting planes. 

Notice that in the linear approximation ()18)) the field gradient, dE^/dx, is independent 
of position x. Thus the tune shift experienced by each particle in the bunch is the same 
(a coherent tune shift). However, the exact result ()19)1 demonstrates that the coherence is 
violated and equation (fT^ allows us to estimate the accuracy of the linear approximation. 

4 Magnetic Images 

In the above, we have used electrostatic images. Magnetic images can be treated in much 
the same way. Let the ferromagnetic boundaries be represented by a pair of infinitely wide 
parallel surfaces at x = ±(7. The magnetic field lines must enter the magnetic pole faces 
perpendicularly. For magnetic image fields we distinguish between DC and AC image fields. 
The DC field penetrates the metallic vacuum chamber and reaches the ferromagnetic poles. 
In case of bunched beams the AC fields are of rather high frequency and we assume that 
they do not penetrate the thick metallic vacuum chamber. 

The DC Fourier component of a bunched beam current is equal to twice the average 
beam current qc(3\B -5^, where B is the the Laslett bunching factor. Thus, 

_ ^0 2qN(3c ^ 4 
By,image,DC[X,X) = — B--A{r],r]), (25) 

where rj = x/g and fj = x/g and the function A is of the form (|2H) . 

The contribution of magnetic AC image field due to eddy currents in vacuum chamber 
walls is similar to electric image fields 

By,image,Aci^,^) = £^^(1 " ^) " ^A(<^,5"), (26) 

where the factor (1 — B) accounts for the subtraction of the DC component. 



The magnetic image fields must be added to the direct magnetic field from the bunch 
to meet the boundary condition of normal components at ferromagnetic surfaces. That is, 
the summary magnetic field between the conducting planes is 

^y,tot\XjX) — Ijy + ijy,image,DC ' tiy,image,AC 



|i + (1 _ B)lAi6, 6) - B-A{7j,fj)}. (27) 

^x n q i 



Att L Ix h g 



5 Summary 



We have derived an approximate expressions for electric (|14|) and magnetic (|16j) self-fields 
produced by a relativistic circular bunch with uniform charge density. They show that at 
distances far from the bunch the electromagnetic field coinsides with the field generated by a 
point-like charged paricle. At the same time, in the near region and beyond the bunch tails, 
the fields coincide with the external self- fields of a continuous beam ©-©• 

We re-analyzed the problem of summing the image fields generated by a bunch of charged 
particles moving with a relativistic velocity between infinitely wide parallel conducting planes. 
The exact solution of the problem represented by the structure function of image fields A 
()2U() depending only of the normalized variables. 
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Appendices 
A Image Fields in Vicinity of a Bunch 

Here we derive the main formula ()2U() . 

Let split the contribution of all image fields (|17j) given in braces into two parts, 

(2/i-xi)-i -(2/i + xi)-i -{ih-X2)-^ +(4/i + X2)-^ 
+ (6/i-xi)"^ -{6h + xiy^ -{8h-X2y^ + {8h + X2)~'^ 
+ (lOh-xi)-^ -{lOh + xiy^ -{l2h-X2)~^ + {l2h + X2y^ + ... (28) 

oo oo 

= ^nfc(xi,/i)-^n„(x2,/i), (29) 

k ra 

where 11^ represents the contribution from the negative charged images and 11^ is the contri- 
bution from the positive charged images. Here and hereinafter, indexes k and m are possess 
odd, A;=l,3,5,..., and even, m,=2,4,6,... values. 

An expansion of denominators of H^ and Hm into a power series of small parameters 
^i = xi/h < 1 and S2 = X2//1 < 1 gives 

^. .. ^ _J: 1 _ 2x1 _2^ J^-^ 

iifcixi,nj 2kh-xi 2kh + xi {2khf-x\ h^^{2kf^' ^ ' 

jr ( ,. ^ _^ 1 _ 2X2 _ 2 ^ 5l"-' 

^^m{X2,n) 2mh-X2 2mh + X2 {2mhy - xj /i^(2m)2"- ^"^^ 
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Now it is evident that the space structure of the image fields between planes is characterized 
by a specific function A{5i,52), we term it the structure function, 



with 



A{5i,S2 



Ef^ + 



k m 



h 



Ai6i,62). 



i{2ky {2kY 



OO (- 



51 



(2m) 



+ .. 



(32) 



(33) 



The structure function A depends only on the normalized variables. 
To proceed further, let us define the following auxiliary quantities 



Mj 



(-) 



M 



(+) 



E 

k 

OO 



1 



(2A;)2i 
1 



E 

m 

OO 



1 



(-1) 



n.+l 






(2m)2i 
1 



^ (2n)2i 

n.=l 
OO ^ 

E 



1 (22i-l-l)^2i 

= w — m — '^^^'i' ^^^^ 

1 22^-1 . ^2i 



(2m)2i ^ (2n)2i 22i (2i)! 



5^ 



2j\, 



(35) 



where i?2j are Bernoulli numbers, B2 = 1/6, B4 = —1/30, Bq = 1/42 etc. By adding and 
subtracting the leftmost parts of (|5H) and (|S^ . we express k and tti numerical series of (|5S]) 
in terms of M- and M- . Therefore, we get from (|33|1 



^ 00 



n=l 



(36) 



or after substituting of (|34() - (|35|) in (|36() . we find the following form of the structure function 
generated by the charged bunch. 



^ 00 



\)b{ 



2n-l r2n-l 



^2n 



22'^(2n)! 



IB 



2n| 



(37) 



Using only the linear terms we recover the part derived by Laslett PP (see equation ()18() 

1 



A(x, x,h) = — ■ ei(2x + x). 



(38) 



An inspection of 1)3 7|) shows that the contributions of negative charged images are enhanced 
by the factor 22" — 1, as compared with the contributions from the positive charged images. 
Equation (|37j) also shows that for x in the bunch center, ^2 = and the contributions from 
the positive charged images are vanish. 

At the final step, it is possible to rewrite the infinite series (|37j) in terms of elementary 
trigonometric functions. To do this, recall the relations between the Bernoulli numbers and 
the trigonometric functions [SI ,10] 



~ (22--l)(2z)2" „ , - 

ztan{z) = }2- (^ V \B2n\, z cot(z) = 1 - ^ 



(2z 



,2n 



n=l 



n=l 



(2n)\ 



IB 



2n\ 



(39) 



After some algebraic manipulations and the use of (|39|1 . we get from (|37|) a new exact and 
compact expression of the structure function 



A('^l,'^2) = 2 [l*™(l^l) + I^°*(l'^2) - Y 



(40) 



Now, if we recall that 6i = (x + x)/h = 5 + 5 and 82 = (x — x)/h = 5 — 5, we obtain 



A{6,5) 



vr 



cos(fJ) 



(41) 



-2 sin(f 5) - sin(f (5) 6-6i 

At a first glance, equation (|ifl|) or (|^T|) is singular at 52 = or (5 = 5, respectively. 
However, as we already discussed right after (|57j) . it is not the case. Starting once again from 
pT]) with 62 = and account (|39|) . we get formally 



A(5,5) = ^tan(p). 



(42) 



Equations (|37() . (|4U() and 1)41(1 were derived assuming (^ < 1 and 5 < 1. Therefore, one 
cast doubts on validity of (|1T|) at (5 ~ 1, near the conducting plane. For that reason in the 
next section we re-expand the series ()28() into a power series of new small parameters. 

B Image Fields in Vicinity of a Conducting Plane 

A similar derivation is used to obtain the field structure near a conducting surface. For 
the case under consideration we have to choose new small parameters for the expansion. 
Each bracket in ()28|) we represent in the form (1 it A)~^ and expand in series, recalling that 
at the plane x ^ h, 

h — xi , , ^ h — X2 

Ai = — -^ < 1, and A2 = — ^ « 1. 



h 
In this way, 

{2kh-xi)-^ = [{2k-l)h]-^ 

{2kh + xi)-^ = [{2k + l)h]-^ 

{2mh-X2Y^ = [{2m-l)hY^ 

{2mh + X2r^ = [{2m + l)h]-^ . ^ 

L Im + 1 

Let us introduce the following auxiliary notations 



1 + 


Ai 1 


2k -l\ 
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Ai 1 




2A; + lJ 


1 + 


A2 


2m -1 


1 


A2 



h 



1 1 



n— 1 A n— 1 



^ (2k - 1)" 

n=l ' 



1 1 



A' 



n-l 



h^(2k + l 



^(2A; + 

n=l ^ 



1 1 



E 



-1)' 



-1 A"-l 



h ^ (2m - 1)" 
h ^ (2m + 1)"' 



L 



n, 



k 
00 



n 



.ly-i 



1 



(2A; - 1)J (2/c + l)i 

I'J' ^ {2k - 1)J ^ (2m - l)i 4^ 

2'J ^^ (2/fc + 1)J -^ (2m + l)i ^ 

K III' n — J_ 

I'J ^^ (2fc - IV -^ (2m - IV ^ (2n - iV 



(2m - ly (2m + l)i 



(2n + l)i 



L 



(+) 
2,i 



00 



{2k - ly ^ (2m - \y ^ (2n 
^ (2fc + 1)J "^ ^ (2m + l)i ~ ^ (2n + 1)J 



(43) 
(44) 
(45) 
(46) 

(47) 
(48) 
(49) 
(50) 
(51) 



By simple manipulations with seriess (jlH|) - (|49j) and (|^ - (|HT|) . it is easy to prove that 



r{-)_i r(-) 



r{+)_r(+) -I 



From l|48j) - (|52j) now easy to find 

4+) +4") = [(-i)-i-i]LW + [(-ir-i + i]Ly, 
4+) -4") = [(-i)-i-i]4+)-[(-i)-i + i]4-j + 2. 



(52) 

(53) 
(54) 



Let introduce the image field structure function A in the way similar to (|32)) and rewrite A 
in terms of notations (|47|1 



A(Ai,A2) 



1 r r / (-i)"~ 

n=l k 

oo 

-Ar'[E( 



1 



(2A;-1)" (2fc + l)" 
(-!)"-! 1 



} 



.(2m-l)'^ (2m + 1) 

^ oo 



(55) 



n=l 



To perform the summation in (|55j) . we have to split the series (|55j) into even, n = 2i, and 
odd, n = 2i — 1 parts and substitute (f^ . (fSl]) in equation H55() 



A(Ai,A2) 



^ oo 

i ^ [(Af + Ar)Lii^, - (Ar 1 - Ar+ML(+) 



n=0 



A2n A2n+1 
l,2(n+l) ^2 ^2 



(56) 



With the help of (011), JSOl) and we find 

r2n+l 



L 



{-) 

l,2n+l 



TT 



22(n+l)(2n)! 



LE. 



■2ra| 



L 



{+) 

l,2(n+l) 



and 



£(Ar+A 



2n+l\ 
2 J 



(22{«+l) _ l)7r2("+l) 

2-[2(n + l)]! 
1 



B. 



2(n+l)l) 



n=0 



I-A2' 

where Bn and £'„ are Bernoulli and Euler numbers, respectively. Thus, 

^2 



(57) 



(58) 



A(Ai,A2) = -|2-- ^ + — (A2-A 

^ ' ^ 2 1 4 1- A2 8 ^ 



+ ^ [(Af + Af )Li72;+i - (Af 

n=l 



2n , a2?in,;'(-) _ /A2n+1 _ ^271+1 \^(+) 



l,2(n+l) 

By using the results obtained in the previous section and the decomposition 



}■ 



(59) 



sec 



(-) = E 



n=0 



\E2n \ 



„2n 



:X 



(60) 



equation (|59|) can be finally expressed in terms of trigonometric functions, 

./. .^ 1 f TT ,vr.. vr ,vr.. vr ,vr.. vr ,'7r., 1 1 ,, 

A(Ai, A2) = -|-sec(-Ai) - -tan(-Ai) + -sec(-A2) + -tan(-A2) - YTa^J- (^1) 

9 



The structure function A(Ai, A2), as written in (|(il|) . looks very different from (|4()|) . However, 
it is not difficult to check that by use of the relations Ai = 1 — 5i and A2 = 1 — 62, equation 
dnU transforms in (gHI) or (|IT|) . 

In this way we ensure that equations (|4U() . (|4H) and (|61() are correct and represent exact 
summation of image fields generated by a charged bunch between infinitely wide conducting 
planes. 
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